arXiv:hep-ex/0211013vl 6 Nov 2002 


SOLAR NEUTRINO OBSERVATIONS AT THE 
SUDBURY NEUTRINO OBSERVATORY 


A.W.P. Poon 

Institute for Nuclear and Particle Astrophysics 
Lawrence Berkeley National Laboratory, Berkeley, CA 94720 

Representing the Sudbury Neutrino Observatory Collaboration* 


*Q.R. Ahmad, R.C. Allen, T.C. Andersen, J.D. Anglin, J.C. Barton, E.W. Beier, M. Bercovitch, J. Bigu, 

S. D. Biller, R.A. Black, I. Blevis, R.J. Boardman, J. Boger, E. Bonvin, M.G. Boulay, M.G. Bowler, 

T. J. Bowles, S.J. Brice, M.C. Browne, T.V. Bullard, G. Biihler, J. Cameron, Y.D. Chan, H.H. Chen, 
M. Chen, X. Chen, B.T. Cleveland, E.T.H. Clifford, J.H.M. Cowan, D.E. Cowen, G.A. Cox, X. Dai, 
E. Dalnoki-Veress, W.E. Davidson, RJ. Doe, G. Doucas, M.R. Dragowsky, C.A. Duba, EA. Duncan, 
M. Dunford, J.A. Dunmore, E.D. Earle, S.R. Elliott, H.C. Evans, G.T. Ewan, J. Earine, H. Eergani, 
A.R Eerraris, R.J. Eord, J.A. Eormaggio, M.M. Eowler, K. Erame, E.D. Erank, W. Erati, N. Gagnon, 

J. V. Germani, S. Gil, K. Graham, D.R. Grant, R.L. Hahn, A.L. Hallin, E.D. Hallman, A.S. Hamer, 
A.A. Hamian, W.B. Handler, R.U. Haq, C.K. Hargrove, RJ. Harvey, R. Hazama, K.M. Heeger, 
W.J. Heintzelman, J. Heise, R.L. Helmer, J.D. Hepburn, H. Heron, J. Hewett, A. Hime, J.G. Hykawy, 
M.C.R. Isaac, R. Jagam, N.A. Jelley, C. Jillings, G. Jonkmans, K. Kazkaz, RT. Keener, J.R. Klein, 

A. B. Knox, R.J. Komar, R. Kouzes, T. Kutter, C.C.M. Kyba, J. Law, I.T. Lawson, M. Lay, H.W. Lee, 

K. T. Lesko, J.R. Leslie, I. Levine, W. Locke, S. Luoma, J. Lyon, S. Majerus, H.B. Mak, J. Maneira, 
J. Manor, A.D. Marino, N. McCauley, D.S. McDonald, A.B. McDonald, K. McEarlane, G. McGregor, 
R. Meijer, C. Mifflin, G.G. Miller, G. Milton, B.A. Moffat, M. Moorhead, C.W. Nally, M.S. Neubauer, 
EM. Newcomer, H.S. Ng, A.J. Noble, E.B. Norman, V.M. Novikov, M. O’Neill, C.E. Okada, R.W. Oiler- 
head, M. Omori, J.L. Orrell, S.M. Oser, A.W.R. Boon, T.J. Radcliffe, A. Roberge, B.C. Robertson, 
R.G.H. Robertson, S.S.E. Rosendahl, J.K. Rowley, V.L. Rusu, E. Saettler, K.K. Schaffer, M.H. Schwen- 
dener, A. Schiilke, H. Seifert, M. Shatkay, J.J. Simpson, C.J. Sims, D. Sinclair, R. Skensved, A.R. Smith, 
M.W.E. Smith, T. Spreitzer, N. Starinsky, T.D. Steiger, R.G. Stokstad, L.C. Stonehill, R.S. Storey, 

B. Sur, R. Tafirout, N. Tagg, N.W. Tanner, R.K. Taplin, M. Thorman, RM. Thornewell, RT. Trent, 
Y.I. Tserkovnyak, R. Van, R.G. Van, C.J. Virtue, C.E. Waltham, J.-X. Wang, D.E. Wark, N. West, 
J.B. Wilhelmy, J.E Wilkerson, J.R. Wilson, R. Wittich, J.M. Wouters, M. Yeh 



ABSTRACT 


The Sudbury Neutrino Observatory (SNO) is a 1000-tonne heavy wa¬ 
ter Cherenkov detector. Its usage of D 2 O as target allows the simultaneous 
measurements of the v^. flux from decay in the Sun and the total flux 
of all active neutrino species through the charged-current and the neutral- 
current interactions on the deuterons. Assuming the standard ®B shape, 
the z/e component of the ®B solar neutrino flux is measured to be 0 e= 
l-76lo;o5(stat.)+Q Qg (syst.) x 10 ® cm“^s“^ for a kinetic energy threshold of 
5 MeV. The non-z/gcomponent is found to be = S.dltg; 45 (stat.)lQ 45 (syst.) x 
10® cm“^s“^. This S.Sa difference provides strong evidence for z/g flavor 
transformation in the solar neutrino sector. The total active neutrino flux is 
measured with the neutral-current reaction at a neutrino energy threshold of 
2.2 MeV. This flux is determined to be 0 nc= S.OOt®; 43 (stat.)lo 4 ® (syst.) x 
10® cm“^s“^, and is consistent with solar model predictions. Assuming an 
undistorted ®B spectrum, the night minus day rate is 14.0±6.3(stat.)ti;4(sys.)% 
of the average rate in the charged-current channel. If the total active neu¬ 
trino flux is constrained to have no asymmetry, the night-day asymmetry in 
the z/e flux is found to be 7.0±4.9(stat.)l}; 2 (sys.)%. A global analysis of all 
the available solar neutrino data in terms of matter-enhanced oscillations of 
two active flavors strongly favors the Large Mixing Angle (LMA) solution. 



1 Introduction 


For more than 30 years, solar neutrino experimentsS’ii’iili have been observing fewer 
neutrinos than what are predicted by the detailed modelsSi of the Sun. This deficit is 
known as the Solar Neutrino Problem. A comparison of the predicted and observed 
solar neutrino fluxes for these experiments are shown in Table These experiments 
probe different parts of the solar neutrino energy spectrum, and show an energy depen¬ 
dence in the observed solar neutrino flux. These observations can be explained if the 
solar models are incomplete or neutrinos undergo flavor transformation while in transit 
to the Earth. 

The Sudbury Neutrino Observatory! was constructed to resolve this solar neutrino 
puzzle. It is capable of making simultaneous measurements of the electron-type neu¬ 
trino (i^e) flux from decay in the Sun and the flux of all active neutrino flavors 
through the following reactions: 


Oe + d 

P + P + e~ 

(CC) 

u^ + d - 

p + n + Ux 

(NC) 

Ux + e~ - 

z^x + e" 

(ES) 


The charged-current (CC) reaction on the deuteron is sensitive exclusively to z/g, and 
the neutral-current (NC) reaction has equal sensitivity to all active neutrino flavors (z/a,; 
a; = e, /r, r). Elastic scattering (ES) on electron is also sensitive to all active flavors, but 
with reduced sensitivity to z/^ and Vr- 

Because of its equal sensitivity to all active neutrinos, the neutral-current measure¬ 
ment can determine the total neutrino flux, hence resolving the Solar Neutrino Problem, 
even if neutrinos oscillated SNO is currently the only experiment that can simultane¬ 
ously observe the disappearance of 1 /^ and the appearance of another active species. 
This is illustrated in Eigure [I]. Another feature of the neutral-current interaction is its 
low kinematic threshold. By efficiently counting the free neutrons in the final state of 
the neutral-current reaction, the total active ®B neutrino flux can be inferred for neutri¬ 
nos with energy above the 2.2-MeV kinematic threshold. 

Recent resultsSS from the measurements of the solar ®B neutrino flux by the SNO 
detector using the CC, NC and ES reactions are presented in this paper. The resultslii 
of a measurement of the day-night asymmetry of the neutrino event rates, which is 
predicted under certain neutrino oscillation scenarios, are also presented. Einally, the 
physics implications of these observations will be discussed. 


Experiment Measured Flux SSM FluxI 

HomestakeS 2.56±0.16(stat.)±0.16(sys.) SNU 7.6tl;?SNU 

SAGskl 70.8iy(stat.)l|^ SNU 128t^ SNU 

GallexB 77.5±6.2(stat.)lt?(sys.) SNU 128t? SNU 

GNoi 65.Sll%^istat)1:liisys.) SNU 1281? SNU 

Kamiokandey 2.80U0.19(stat.)±0.33(sys.)x 10® 5.05 xlO® 

Super-Kamiokandei 2.35±0.03(stat.)lg:g^(sys.)x 10^ 5.05x10^ (lloae) 


Table 1: Summary of solar neutrino observations at different solar neutrino deteetors. 
The measured fluxes at the radioehemieal experiments are measured in Solar Neutrino 
Unit (SNU), whieh is defined as 1 eapture per seeond per 10^® target atoms. The 
Kamiokande and Super-Kamiokande measurements are in units of 10® em“^ s“l 
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Figure 1: Using the measured solar neutrino fluxes from different reaetion ehannels to 
provide “smoking gun” evidenee of neutrino flavor transformation. 














2 The Sudbury Neutrino Observatory 


2.1 Physical Description of the SNO Detector 

SNoi is an imaging water Cherenkov detector located in the Creighton #9 mine, owned 
by the International Nickel Company (INCO) near Sudbury, ON, Canada. A barrel¬ 
shaped cavity with a height of 34 m and a diameter of 22 m was excavated at a depth 
of 2092 m (or 6000 meters of water equivalent) to house the detector. The muon flux 
traversing the detector is 67 day“^. 

Figure ^ shows a cross-sectional view of the SNO detector. The neutrino detec¬ 
tion medium is 1000 metric tons of 99.92% isotopically pure D 2 O contained in a 12-m 
diameter acrylic sphere. The acrylic vessel is constructed out of 122 ultraviolet trans¬ 
mitting acrylic panels. This sphere is surrounded by 7000 metric tons of ultra-pure H 2 O 
contained in the cavity. This volume of H 2 O shields the detector from high energy 7 
rays and neutrons originating from the cavity wall. A 17.8-m diameter stainless steel 
structure supports 9456 20-cm diameter inward-facing photomultiplier tubes (PMTs). 
A non-imaging light concentrator is mounted on each PMT, extending the total pho¬ 
tocathode coverage to 55%. An additional 91 PMTs are mounted facing outward on 
the support structure to serve as a cosmic-ray veto. To cancel the vertical components 
of the terrestrial magnetic field, 14 horizontal magnetic compensation coils were built 
into the cavity wall. The maximum residual field at the PMT array is <19/rT, and the 
reduction in photo-detection efficiency is about 2.5% from the zero-field value. 

A physics event trigger is generated in the detector when there are 18 or more PMTs 
exceeding a threshold of ~0.25 photo-electrons within a coincidence time window of 
93 ns. All the PMT hits registered in the ~420 ns window after the start of the coinci¬ 
dence time window are recorded in the data stream. This widened time window allows 
scattered and reflected Cherenkov photons to be included in the event. The mean noise 
rate of the PMTs is ~500 Hz, which results in ~2 noise PMT hits in this 420 ns win¬ 
dow. The instantaneous trigger rate is about 13-20 Hz, of which 5-8 Hz are physics 
triggers. The remaining triggers are diagnostic triggers for monitoring the well being 
of the detector. The trigger efficiency reaches 100% when the PMT multiplicity (N/jjt^) 
in the event window is >23. 



Figure 2: A cross-sectional view of the SNO detector. The outer geodesic structure is 
the PMT support (“PSUP”). 












































2.2 Solar Neutrino Physics Program at SNO 

The solar neutrino physios program at SNO is designed to exploit its unique NC oapa- 
bility. Beoause the result of this NC measurement is a definitive statement on the Solar 
Neutrino Problem and flavor transformation of solar neutrinos, the SNO experiment 
has implemented a plan to make three independent NC measurements of the total ®B 
aetive neutrino flux. 

The first NC measurement was made with a pure D 2 O target. The free neutron 
from the NC interaetion is thermalized, and for ~30% of the time, a 6.25-MeV 7 ray 
is emitted following the neutron eapture by a deuteron. Only ~50% of the 6.25-MeV 
observed photopeak is above the neutrino analysis threshold, yielding a deteetion ef- 
fieieney of ~15%. The results from this NC measurement is presented in this paper. 
The seeond NC measurement is being made with 2 tonnes of NaCl added to the D 2 O. 
The free neutron is readily eaptured by ^®C1 in this deteetor eonfiguration, and a eas- 
eade of 7 rays with a total energy of 8.6 MeV follow. The neutron deteetion effieieney 
is signifieantly enhaneed, and ~45% of the NC events have a deteetable signal above 
the analysis threshold. This phase of the experiment is seheduled to eomplete by the 
end of 2002. In the third NC measurement, diserete ^He proportional eounters will be 
installed inside the D 2 O volume.0 The neutron deteetion effieieney of the proportional 
eounter array is ~40%. In this deteetor eonfiguration, the deteetion of the CC and the 
NC signals are deeoupled, and the eovarianee of the CC and NC signals that appear in 
the first two deteetor eonfigurations is eliminated in this ease. 

3 Solar Neutrino Analysis at SNO 

The data presented in this paper were reeorded between November 2,1999 and May 28, 
2001. The eorresponding livetime is 306.4 days. During this data aequisition period, 
the Sun was above the deteetor’s horizon (“day” data set) for 128.5 days, and below the 
deteetor horizon (“night” data set) for 177.9 days. The target was pure D 2 O during this 
period. Figure ^ summarizes the analysis proeedure for extraeting the CC, NC and the 
ES event rates in the SNO deteetor. In the following, eaeh step in the analysis flow is 
diseussed in detail. 


Data 



Figure 3: Simplified flow ehart of solar neutrino analysis at SNO. 

3.1 Data Reduction and Data Loss 

After removing all the deteetor diagnostie triggers from the data stream, a signifieant 
portion of the remaining events are instrumental baekgrounds. Examples of these baek- 
grounds inelude eleetrieal diseharges in the PMTs (“flashers”) and light emission from 
diseharges in insulating deteetor materials. Data reduetion sehemes were developed to 
remove these baekgrounds. 

3.1.1 Pass 0 Cut 

The instrumental baekgrounds have eharaeteristie PMT time and eharge distributions 
that are signifieantly different from Cherenkov light, and ean be eliminated using euts 
based on these distributions. For example, the diseharge light emitted from a flasher 
PMT is detected across the detector ~70 ns after the initial discharge is registered. 
Some of these light-emitting instrumental backgrounds are localized near the water 
piping near the top of the detector. Veto PMTs were installed in this region in order 
to enhance the rejection efficiency of these non-Cherenkov events. Interference in the 








electronics system can produce false events. Most of the registered electronic chan¬ 
nel charges in these interference events are near the pedestal, and can be removed by 
a cut on the mean charge of the fired PMTs. Some of these electrical discharge or 
electronic interference background events also have different event-to-event time cor¬ 
relations from physics events, and time correlation cuts are used to remove these events. 
Two independent instrumental background rejection schemes are used. An event-by- 
event comparison of the data sets reduced by these two schemes shows a difference of 
< 0 . 2 %. 

The physics loss due to these instrumental background cuts is calibrated with a trig¬ 
gered 6.13-MeV 7 -ray sourcelli and a triggered ®Li 13-MeV endpoint [3 source^ 
deployed to the D 2 O and H 2 O volumes. Further tests of the 13ihits dependence in the 
cuts are performed with an isotropic light source at various intensities. The physics 
acceptance of the instrumental background cuts, weighted over the fiducial volume, is 
measured to be 0.9966lo;ooo2. Instrumental background rejection is well over 99% at 
this stage. 

3.1.2 High Level Cuts 

After passing the instrumental background cuts, all events with N/ij4s>30 (~3.5 MeV 
electron energy) are reconstructed. Once the event reconstruction information becomes 
available after the reconstruction, several high level physics cuts are applied to the 
Pass 0-reduced data set to further reduce the instrumental backgrounds. The efficiency 
of the reconstruction algorithm is calibrated with the and the ®Li sources, and is 
found to be 0.9985±0.0015 for neutrino events originated from the central 550 cm of 
the detector. 

The high level cuts test the hypothesis that each event has the properties of elec¬ 
tron Cherenkov light. The reconstruction figure-of-merit cuts test for the consistency 
between the time and angular expectations for an event fitted to the location of the re¬ 
constructed vertex and that based on the properties of Cherenkov light and the detector 
response. 

Two parameters are used to further characterize Cherenkov light. The average open¬ 
ing angle between two hit PMTs ((%)), measured from the reconstructed vertex, is used 
to determine whether the topology of an event is consistent with Cherenkov light. The 
in-time ratio (ITR) is the ratio of the number of hit PMTs within an asymmetric time 
window around the prompt light peak to the number of calibrated PMTs in the event. 


Channel 

Pass-0 

Reconstruction 

High Eevel Cuts 

Total 

CC 

O.yyOO 0.0002 

0.9985±0.0015 

0.9906T0.0005 

0 . 9857 ^°^ 

ES 

U.yyOO 0.0002 

0.9985±0.0015 

0.9903±0.0007 

o.9854!!:K 

NC 

O.yVOO 0.0002 

0.9985T0.0015 

0.9821T0.0011 

0.9772^0^ 


Table 2: Physics acceptance at each of the instrumental background removal stages. A 
kinetic threshold of 5 MeV and a fiducial volume of the inner 550 cm of the detector 
are assumed. 

Figure ^ shows the correlations between 9ij and ITR for instrumental backgrounds and 
Cherenkov light events. As shown in the figure, this two dimensional cut has very high 
instrumental background rejection efficiency. 

The total signal loss from the Pass 0 and the high level cuts are calibrated with the 
and the ®Li sources. Because of the difference in the energy spectrum and spatial 
distribution for CC, ES and NC events, there is a slight difference in the data loss due to 
the Cherenkov likelihood cut for the different detection channels. Table ^summarizes 
the physics acceptance after the Pass-0, reconstruction and high level cuts for the three 
neutrino detection channels. 

The residual instrumental background contamination in the neutrino signal after the 
series of instrument background cuts is estimated by a bifurcated analysis, in which 
the signal contamination is obtained from cross calibrating the background leakage of 
two groups of orthogonal cuts. For the same fiducial volume and energy threshold, the 
instrumental background contamination is estimated to be <3 events (95% C.L.), or 
0.1% of the final neutrino candidate data set. Table ^ summarizes the sequence of cuts 
that are used to reduce the raw data to 2928 neutrino candidate events. 

3.2 Reconstruction and Energy Calibration 

After all the instrumental background cuts have been applied to the data, energy and 
fiducial volume cuts are employed to reduce the physics backgrounds in the neutrino 
candidate data set. These physics backgrounds include low energy backgrounds from 
radioactive decays in the natural and ^^^Th chains, and other backgrounds induced 
by cosmic rays. Because of progressively lower radioactivity towards the central region 
of the detector, one can place a fiducial volume cut and an energy threshold sufficiently 
high that most of the low energy backgrounds are removed from the neutrino candidate 
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Figure 4: Separation of instrumental backgrounds and Cherenkov light events using the 
high level cuts. 


Analysis step 

Number of events 

Total event triggers 

450 188 649 

Neutrino data triggers 

191 312 560 

PMT hit multiplicity (Nhus) >30 cut 

10 088 842 

Instrumental background (Pass 0) cuts 

7 805 238 

Cherenkov likelihood cuts 

3 418 439 

Fiducial volume cut 

67 343 

Energy threshold cut 

3 440 

Cosmic-induced background subtraction 

2 928 

Total events 

2 928 


Table 3: Number of candidate events remained after each data reduction step 







data set. These euts rely upon event reeonstruetion and deteetor energy ealibration, 
whieh we will diseuss in the following. 

3.2.1 Reconstruction 

The ealibrated times and positions of the fired PMTs are used to reeonstruet the ver¬ 
tex position and the direetion of the partiele. Two different reeonstruetion algorithms 
were developed. An event-by-event eomparison shows exeellent agreement between 
the data sets reeonstrueted by these two algorithms. The neutrino data presented in 
this paper are reeonstrueted using a maximum likelihood teehnique whieh uses both 
the time and angular eharaeteristies of Cherenkov light. Vertex reeonstruetion aeeu- 
raey and resolution for eleetrons are measured using Compton eleetrons from the 
7 -ray souree, and their energy dependenee is verified by the ^Li f3 souree. Compton 
seattered eleetrons from a 6.13-MeV 7 ray are preferentially seattered in the forward 
direetion relative to the ineident 7 -ray direetion. Thus they provide information about 
the aeeuraey of the direetion reeonstruetion. In order to minimize the effeet of finite 
vertex resolution on this angular resolution measurement, only events that are re¬ 
eonstrueted to more than 150 em from the souree are used in the measurement. At the 
energy (~5.5 MeV total eleetron energy), the vertex reeonstruetion resolution is 
16 em and the angular resolution is 26.7°. Reeonstruetion-related systematie uneertain- 
ties to the solar neutrino flux measurement are 1^8% for the CC ehannel and for 
the NC ehannel. 

3.2.2 Energy Estimator 

Figure ^ shows the relationship between the deteetor ealibration program and event-by¬ 
event energy estimation in the analysis. Optieal ealibration is obtained using a near- 
isotropie souree of pulsed laser light0’@ at 337, 365, 386, 420, 500 and 620 nm. The 
light souree is deployed to loeations aeeessible by the souree manipulator system on 
two orthogonal planes in the D 2 O, and on a linear grid in the H 2 O. Optieal parameters 
of different optieal media in the deteetor are obtained at these wavelengths .0 The 
attenuation lengths in D 2 O and H 2 O are found to be near the Rayleigh seattering limit. 
These optieal parameters are inputs to the Monte Carlo/energy estimator engine. 

The triggered souree is used to provide the absolute energy ealibration. The 
deteetor energy response to the photopeak of this souree provides a normalization to 
the PMT photon eolleetion effleieney used in the Monte Carlo, and establish the abso- 
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Figure 5: Calibration of the SNO detector and event-by-event energy estimator. 


lute energy calibration. A long-term stability study of the detector response to the 
source shows a linear drift of 1.3% year“^. The cause of this effect is under investiga¬ 
tion, and a drift correction is applied to the event-by-event energy estimator. 

This tuned Monte Carlo is then used to make predictions for the energy response 
to different calibration sources. The pT source generates 19.8-MeV 7 rays through 
the ^H(p, 7 )^He reaction]^ and is used to check the linearity of the energy response 
beyond the endpoint of the neutrino energy spectrum. To probe the regions that 
are inaccessible to the calibration source deployment system, the ^®^Cf fission neutron 
source provides an extended distribution of 6.25-MeV 7 rays from d(n, 7 )t. Figure ^ 
shows a comparison of the Monte Carlo predictions and the detector responses to these 
sources. 

The energy estimator uses the same input parameters (e.g. optical parameters) as 
the Monte Carlo. It assigns an effective kinetic energy Teg to each event based upon its 
position, direction and the number of hit PMTs within the prompt (unscattered) photon 
peak. For an electron of total energy the derived detector energy response can be 
parameterized by a Gaussian: 


Ee) 


1 

,— -exp 

V^aE{Ee) 


1 f Ees-E^ y 
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Figure 6 : Comparison of the Monte Carlo predicted responses to different calibrated 
sources. 

where i?eff=^efr+me, and the energy resolution is given by 

aE{Ee) = -0.0684 + 0.33i/e^^^ + 0.0425(Ee - me) MeV. 

The systematic uncertainty on this absolute energy calibration is found to be ±1.2%, 
which results in neutrino flux uncertainties of t|; 2 % for the CC channel and for 
the NC channel. This is the most significant systematic uncertainty in the flux mea¬ 
surement. Other energy related systematic uncertainties to the flux include the energy 
resolution and the energy scale linearity. A summary of the systematic uncertainties 
can be found in Table 0. 

A second energy estimator using Nhus is employed for validation purposes. These 
two energy estimators give consistent results in the neutrino flux measurement. 

3.3 Neutron Efficiency 

For the NC measurement of the solar neutrino flux, the probability that a neutron pro¬ 
duced in the NC interaction will capture on deuterium and the detection efficiency of 
the 6.25-MeV 7 ray emitted following the neutron capture have to be determined. The 
















capture probability depends on the isotopie abundances of various nuelei in the heavy 
water and the relevant eapture eross seetions for eaeh. Near the boundary of the AV, 
the eapture cross sections for neutrons in acrylie and light water beeome important as 
well. The deteetion effieieney of the 6.25-MeV 7 ray ean be readily dedueed from the 

ealibration data, sinee the latter produees predominantly 6.13-MeV 7 rays. 

Three different methods are used to calibrate the eapture effieieney: a method whieh 
direetly eounts the number of neutrons measured with a ^^^Cf souree and eompares to 
souree expeetations, a purely analytie ealculation eompared to Monte Carlo, and a mul- 
tiplieity measurement which extracts the effieieney statistically from the distribution of 
the number of deteeted neutrons per fission. 

The primary approaeh to measuring the eapture effieieney is the “direct counting” 
method, in whieh the number of neutrons deteeted during a souree run is eompared 
to the total number expeeted to be generated by the source based on the known decay 
rate. This analysis was performed on ^^^Cf source data taken at various loeations in the 
deteetor. Beeause the ^®^Cf souree is not a triggered souree, non-neutron baekgrounds 
associated with the source and the deteetor were handled using reeonstruction and in¬ 
strumental baekground euts, with a proper aeeounting of the assoeiated aceeptanee of 
these euts in the effieieney measurement. In Figure ^ the neutron eapture effieieney 
on deuterium as a function of radial distance from the eenter of the deteetor is shown. 
The eapture effieieney on deuterium for NC events is found to be 29.90 ± 1.10%, and 
when the fidueial volume eut of 550 em and energy threshold cut of >5.0 MeV are 
ineluded, the overall detection efficiency is 14.38 ± 0.53%. The uncertainties on these 
numbers eontribute a relative uneertainty of 3.68% in the extraetion of the NC signal in 
the present analysis. 

The analytie ealeulation uses known eross seetions and isotopic abundances of the 
D 2 O to ealeulate the expeeted effieieney, and propagates uneertainties on the miero- 
scopic parameters through the entire ealeulation. The largest uneertainties in the ana¬ 
lytic model are on the fraetion of H 2 O in the heavy water, the eapture eross seetions 
on deuterium and ^® 0 , and the diffusion parameter which governs the total number of 
neutrons eseaping the D 2 O volume. For the NC signal, the total uneertainty determined 
in this way is 3%. 

The third approach is provided primarily as a eheck beeause of its limited sensi¬ 
tivity. In this approaeh the multiplieity distribution of deteeted neutrons is fit with the 
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Figure 7: Final capture efficiencies measured with the direct counting method as a 
function of radial position, and the residuals from a fit to the radial dependence. 

















Component 

232Th 

238^ 


(g/g) 

(g/g) 

D 2 O 

3.7x10"^^ 

4.5x10“^^ 

H 2 O 

3.7x10-^^ 

4.5x10-13 

AV 

1.9x 10-^2 

3.6x10-12 


Table 4: The target radio-purity levels for different eomponents in the SNO deteetor. 
expression 

_ ^ ^ exp A^exp"^ 

2N'Ka'^ N\ 

where P{d) is the probability to deteet d neutrons per unit time given that r neutrons 
were generated aeeording to Poisson statisties for the ^®^Cf fission rate A and the known 
^®^Cf multiplieity per fission (mean /r = 3.79 ± 0.006 and width cx^ = 1.57 ± 0.02). 
In this ease, both the souree fission rate A and the neutron deteetion effieieney e are 
derived from a maximum likelihood fit. The results from this approaeh and the direet 
eounting approaeh are eonsistent. 

3.4 Low Energy Backgrounds 

Low levels of ^^^Th and ean be found naturally in all materials. The SNO deteetor 
has been designed with very stringent radiopurity targets for different eomponents in the 
deteetor. Table ^ lists the ^^^Th and target levels for the D 2 O, aerylie vessel (AV), 
and H 2 O. At these radiopurity levels, the baekground to the NC signal is approximately 
10% of the NC signal predieted by the Standard Solar Model. 

Radioaetive deeays of the daughters in the natural ^^^Th and ehains are the 
dominant baekgrounds in the neutrino signal window. These deeays ean eontribute to 
the baekgrounds in two different ways. 

A 7 ray with energy in exeess of the binding energy of the deuteron (=2.2 MeV) 
ean potentially photodisintegrate a deuteron: 

7 + d —n + p 

The free neutron in the final state is indistinguishable from that in the NC signature. In 
the natural ^^^Th and the ehains, 7 rays that ean photodisintegrate the deuterons are 
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Figure 8: Low Energy Baekground Analysis. The teehniques that are used in under¬ 
standing these baekgrounds are shown. 

emitted in the deeays of To measure the eontribution of this photodis¬ 

integration baekground to the NC signal, it is neeessary to determine the eoneentration 
of the different isotopes in the deteetor. “Ex-situ” radioehemieal assays and ”in-situ” 
Cherenkov light pattern reeognition teehniques are used to measure the eontribution 
from this photodisintegration baekground. 

A small fraetion of deeays with large Q values (e.g. 208 fj-om the ^^^Th ehain and 
^^"^Bi from the ehain) inside the fidueial volume might generate enough light in 
the deteetor and register a ealibrated energy above the analysis threshold. Deeays that 
are originated from outside the fidueial volume might get reeonstrueted to within the 
fidueial volume due to finite reeonstruetion resolution. Eow energy radioaetive sourees 
and simulations are used to understand this elass of Cherenkov tail baekgrounds. 

Eigure shows the general strategy of the low energy baekground analysis. In the 
following, we shall diseuss the analysis teehniques mentioned above in more details. 

3.4.1 Photodisintegration Background 


D 2 O and H 2 O Radioactivity 




Several techniques were developed to measure the ^^^Th and concentration in 
the D 2 O and the H 2 O. These techniques sample different daughters in the ^^^Th and 
chains, and can be broadly classified into two categories: 

• Ex-situ: These are techniques that involve chemical removal of radioisotopes 
from the D 2 O and the H 2 O, and their decays are counted in a system external to 
the SNO detector. These techniques include the extraction of Ra isotopes using 
MnOx beads (the “MnOx” techniqueil), and the extraction of Ra, Th and Pb iso¬ 
topes using HTiO membranes (the “HTiO” technique^). Typical assays circulated 
approximately 400 tonnes of water through the extraction media. These tech¬ 
niques provide isotopic identification of the decay daughters and contamination 
levels in the assayed water volumes, presented in Fig. |T^(a). Secular equilibrium 
in the decay chain is broken by the ingress of long-lived (3.8 day half-life) 
^^^Rn in the experiment. Measurements of this background are made by period¬ 
ically extracting and cryogenically concentrating ^^^Rn from the chain (the 
“Rn assay”@). The Rn results are presented (as mass fractions in g(U)/g(D20)) in 
Fig. nib). 

The MnOx Technique 

In the MnOx technique, D 2 O or H 2 O is passed through polypropylene columns 
that contain beads coated with a manganese oxide compound (MnOx), which ex¬ 
tracts Ra from the flowing water. After a large volume of water has passed through 
the columns, they are removed and dried. The dried column is then attached to a 
gas flow loop on an electrostatic counter (ESC). The Rn produced from Ra decay 
is swept from the columns into the ESC where it decays. The charged Po ions 
from the decay of Rn are carried by the electric field onto an a counter where the 
decays of the Po are detected, and their a energy spectra are collected. For the 
^^^Th chain, the relevant Po a decays are ^^®Po (6.8 MeV a) and ^^^Po (8.8 MeV 
a), whereas the relevant ones for the U chain are ^^®Po (6.0 MeV a) and ^^^Po 
(7.7 MeV a). 

The HTiO Technique 

In this technique, D 2 O or H 2 O is passed through hydrous titanium oxide (HTiO) 
trapped on filtration fibers. The HTiO ion-exchanger is first deposited onto a 


microfiltration membrane. Then the eolumns eontaining the loaded filters are used 
to extraet ^^^Ra from the Th ehain and ^^®Ra from the U ehain from a large volume 
of D 2 O or H 2 O. After extraetion, Ra is eluted with nitrie aeid, and subsequently 
eoneentrated down to ~10 ml of eluate. This is then mixed with liquid seintillator 
and eounted using (3 — a delayed eoineidenee eounters.0 For the ^^^Th ehain, 
the eoineidenees of the /3-deeay of ^^^Bi and the a-deeay of ^^^Po are eounted, 
whereas the eoineidenees of the /3-deeay of and the a-deeay of ^^“^Po are 
eounted for the ehain. 

Rn Assays 

Measurements of ^^®Ra eoneentration in the D 2 O and the H 2 O in the MnOx and 
the HTiO teehniques are not suffieient to determine the total radioaetive baek- 
ground from the ehain. This is beeause even a small leak of the underground 
laboratory air (~3 pCi// of ^^^Rn) ean lead to signifieant disequilibrium between 
^^®Ra and ^^"^Bi. The Rn assay teehnique was developed to taekle this problem. In 
this teehnique, water drawn from diserete sample points in the deteetor is flowed 
through a degasser to liberate Rn. The Rn is purified and eolleeted in a eryo- 
genie eolleetor. The subsequent a deeays are eounted in a Lueas eell seintillator 
(ZnS) ehamber on a 2.54-om diameter photomultiplier tube. Sinee there is a de¬ 
lay of many ^^°Rn lifetimes between the preparation of the Lueas eells and their 
subsequent eounting, this method is sensitive only to ^^^Rn deeays. 

• In-situ: This is a teehnique that uses pattern reeognition on the Cherenkov light 
distribution to determine the eoneentration of ^^^Th and in the D 2 O and the 
H20.@@ The 208 -pi decay (from the ^^^Th ehain) has a Q value of ~4.9 MeV, 
and the ^^“^Bi deeay (from the ehain) has a Q value of 3.27 MeV. Almost 
every deeay emits a 2.614 MeV 7 , several low energy 7 ’s and a (3 with an 
endpoint of ~ 1-1.8 MeV, whereas there is a unique braneh in the ^^"^Bi deeay 
that produees a single (3 with an endpoint energy of 3.27 MeV. A single parti- 
ele produees more Cherenkov photons than multiple partieles with the same total 
energy beeause of the Cherenkov kinetie threshold. Also the light isotropy is dif¬ 
ferent for these two oases beeause multiple partieles generate a more isotropio 
light pattern (multiple Cherenkov eones), whereas a single partiele gives a more 
direotional light distribution due to light emission within a single Cherenkov eone. 
Therefore, by seleoting a set of energy and fiduoial volume outs, it is possible to 


separate decays statistically by their differences in light isotropy. 

The light isotropy parameter that has been developed ({Oij)) is the average an¬ 
gle between all hit PMT pairs within the prompt light time window in an event. 
Figure ^ shows the difference in 9ij between ^os decays. In this in- 

situ analysis, Cherenkov events fitted within 450 cm from the detector center and 
extracted from the neutrino data set provide a time-integrated measure of these 
backgrounds over the same time period of the neutrino analysis. Statistical sepa¬ 
ration of the 208 and ^^"^Bi events is obtained by a maximum likelihood fit of the 
Cherenkov signal isotropy. 

Results from the ex situ and in situ methods are consistent with each other as shown 
on the right hand side of Figs. [^a) and pDKb). For the ^^^Th chain, the weighted 
mean (including additional sampling systematic uncertainty) of the two determinations 
is used for the background analysis. The chain activity is dominated by Rn ingress 
which is highly time dependent. Therefore the in situ determination was used for this 
activity as it provides the appropriate time weighting. For the present data set, the time- 
averaged equivalent equilibrium and ^^^Th concentration in the D 2 O are found to 
be 


^^^Th ; 1.63 ± 0.58 X IQ-^^gTh/gDaO 
238U : 17.8^11 X 10-^® gU/gDaO. 

Once the 232xh and concentrations are measured, the photodisintegration back¬ 
ground can be determined through Monte Carlo calculations that use the photodisin¬ 
tegration cross section, properties of 7 -ray and neutron transport, and the detector re¬ 
sponse function as inputs. The expected number of observed photodisintegration neu¬ 
trons are 18.4±6.5 counts for the ^^^Th background and 25.9t|;3 counts for the 
background in the pure D 2 O phase of the experiment. 

Daughters in the ^^^Th and chain in the H 2 O and the acrylic vessel (AV) can 
still produce a photodisintegration background. This happens when the 7 ’s with an en¬ 
ergy >2.2 MeV enter the D 2 O volume and photodisintegrate a deuteron. Ex situ assays 
and in situ pattern recognition techniques described above were used to measure the 
^^^Th and concentrations in the H 2 O. The photodisintegration background inside 
the fiducial volume due to radioactivity in the H 2 O was then determined. The time- 
averaged equivalent equilibrium and ^^^Th concentration in the H 2 O are found to 



Figure 9: Difference in 6ij between and ^^^Bi decays. 


be 


: 14.2 ± 6.6 X IQ-^^gTh/gHsO 
: 75.5 ± 33.0 X 10"^^ gU/gHaO. 

At these concentrations, the number of observed photodisintegration neutrons in the 
fiducial volume are determined to be 5 . 6 ^ 2!2 counts for the ^^^Th chain and S.blaig counts 
for the chain. 

Acrylic Vessel Radioactivity 

Prior to the construction of the acrylic vessel, neutron activation analyses were done 
on samples of the acrylic panels used in the vessel. These results, along with that of a 
dust density measurement on the vessel subsequent to its construction, indicate that the 
radioactive load of the acrylic vessel to be 7.5±1J 3 fig of ^^^Th and 15±15 fig of 
Studies of the Cherenkov data revealed a radioactive hotspot on the acrylic vessel. The 
radioactivity of this hotspot is equivalent to fig of ^^^Th. Combining the results 

from radioassays and the hotspot analysis, the expected number of observed photodisin¬ 
tegration background events due to radioactivity in the acrylic vessel is 15.8lg § counts. 





Figure 10: Thorium (a) and uranium (b) backgrounds (equivalent equilibrium concen¬ 
trations) in the D 2 O deduced by in situ and ex situ techniques. The MnOx and HTiO 
radiochemical assay results, the Rn assay results, and the in situ Cherenkov signal de¬ 
termination of the backgrounds are presented for the period of this analysis on the left- 
hand side of frames (a) and (b). The right-hand side shows time-integrated averages 
including an additional sampling systematic uncertainty for the ex situ measurement. 

































































232Th 

(eounts) 

238^ 

(eounts) 

Total 

D 2 O photodisintegration 

18.4T6.5 

25.911° 

44.31®;? 

AV photodisintegration 

14.2^1® 

1.6T1.6 

15.811° 

H 2 O photodisintegration 


5.61^9 


Atmospherie Es and 
sub-Cherenkov threshold /i’s 

Not applieable 

4 ± 1 

Spontaneous Fission 

Not applieable 

< 1 

^H(q;, a)pn 

Not applieable 

2 ±0.4 

^^0(a,n) 

Not applieable 

< 1 

Terrestrial and reaetor z/’s 

Not applieable 


External neutrons 

Not applieable 

< 1 

Total 


78 ± 12 


Table 5: Summary of neutron baekground {E >5.5 MeV) in the fidueial volume 
{R <550 em.) for the pure D 2 O phase of the experiment 

Table ^ is a summary of the photodisintegration baekground for the data set used in 
this analysis. In this table, we have also listed the eontributions from other sourees of 
neutron baekgrounds. Some of these are eosmie ray events and atmospherie neutrinos. 
To reduee these baekgrounds, an additional neutron baekground out imposed a 250- 
ms veto (in software) following every event in whioh the total number of PMTs whioh 
registered a hit was greater than 60. The oontribution from these additional sourees are 
very small oompared to the photodisintegration baekground. 

3.4.2 Cherenkov Tail Background 

Low energy baekgrounds from Cherenkov events in the neutrino signal region were 
evaluated by using aerylie eneapsulated sourees of U and Th deployed throughout the 
deteetor volume and by Monte Carlo ealeulations. In the following, we diseuss how 
this elass of baekground is determined for different origins of radioaetivity. 

Cherenkov Tail Background from D 2 O Radioactivity 


The Cherenkov tail baekground arising from internal radioaetivity in the D 2 O were 
















determined by a Monte Carlo study. A Monte Carlo study is justified because the de¬ 
tector response in the D 2 O is well calibrated, and the internal radioactivity is well mea¬ 
sured by the ex situ and in situ techniques. In this Monte Carlo calculation, the number 
of Cherenkov background events is normalized to the number of photodisintegration 
(pd) neutrons in the neutrino signal window and are found to be: 

208ti: 0.162 1o!o3o Cherenkov tail event/pd neutron event 

^^■^Bi: 0.670 lo!i 25 Cherenkov tail event/pd neutron event. 

Given the number of detected photodisintegration neutrons (Th: 18.4±6.5, U: 25.9^|3), 
the expected Cherenkov tail backgrounds are: 

232Th: 3.0 1?;!] counts 

17.4 ^^^ 3 ^ counts 
Total: 20.4 counts. 

As a consistency check, this Monte Carlo technique was used to make predictions 
on the energy spectrum of a variety of background sources, including a stainless-steel 
encapsulated Th source and a Rn spike in the detector during a water circulation pump 
failure. Figure [H] compares the predicted energy spectrum for the stainless-steel en¬ 
capsulated Th source run to the actual calibration source data. 

Cherenkov Tail Background from External Radioactivity 

In principle, one can perform Monte Carlo calculations to predict the Cherenkov tail 
background in the fiducial volume for radioactivity from the acrylic vessel, the H 2 O and 
the PMT support geodesic. However, most of the events in the neutrino signal window 
with origins in the regions external to the fiducial volume come from the tail of the 
energy response and the reconstruction response functions. These response functions 
are not as well understood outside the fiducial volume. 

A technique was developed to determine this external (acrylic vessel, H 2 O and PMT 
support geodesic) Cherenkov tail backgrounds in the neutrino signal window using the 
low energy background radial distribution (R) calibrated by low energy acrylic en¬ 
capsulated sources. The principle of this technique is to use the acrylic source data 
to generate the probability distribution functions (PDFs) for the low energy back¬ 
grounds in the AV, H 2 O and the PMT regions. The R^ distribution in these regions in 
the neutrino data was then fitted to a linear combination of these PDFs by the extended 
maximum likelihood method. One difficulty of this technique is the lack of sufficient 


Comparison of Predicted Th Source Spectrum 



Figure 11: Comparison of Monte Carlo predietion with systematie uneertainties to a 
stainless-steel eneapsulated Th souree at the eenter of the D 2 O. The /?’s in the deeays 
are bloeked by the stainless steel shroud, while the 7 ’s emitted in the deeays eseaped 
the stainless steel eontainer. The predietions, shown here as a Tier band, have been 
normalized by the livetime and the souree strength. 












Figure 12: External Cherenkov tail fit at Tes>4 MeV. The extended maximum like¬ 
lihood method was used in the fit, and the band represents the systematic uncertainties. 

statistics at the intended neutrino analysis threshold (Teff=5 MeV). In order to overcome 
this, the analysis was proceeded at a lower energy (Teff=4 MeV) and then extrapolated 
to 5 MeV. The basic assumptions here are that there is no correlation between and 
energy, and that the reconstruction does not get worse with higher energy. The R^ fit 
was performed for the range 1.02 < (i?/600)^ < 2.31. Figure shows the results 
of this R^ fit. The band shown in this plot is the range of the systematic uncertain¬ 
ties. The data points are not centered in this band because the systematic uncertainties 
are not normally distributed. Table ^ is a summary of the Cherenkov tail background 
(Teff>5 MeV) in the fiducial volume (R <550 cm.) for the pure D 2 O phase of the 
experiment. 

Several consistency checks of this external Cherenkov tail analysis were performed. 
The radial PDF for H 2 O background was compared to the radial distribution of a radon 
spike during a water circulation pump failure. This is shown in Figure [1^. The expected 















232Th 

(counts) 

238^ 

(counts) 

Total 

D 20 


17.411^3^ 

lOAhlY 

AV 

Not applicable 

O-J-6.3 

H 2 O 

Not applicable 

9 0+3.9 
^• 0 - 2.8 

PMT 

Not applicable 

16.01JV 

Total 

Not applicable 

45.5t\li 


Table 6 : Summary of Cherenkov tail baekground {E >5.5 MeV) in the fidueial vol¬ 
ume {R <550 em.) for the pure D 2 O phase of the experiment. The high energy 7 
baekground eontribution in the fidueial volume is ineluded in the PMT /Jy entry. 

baekground rate in Table ^ was also eheeked against the Monte Carlo ealeulations and 
the results of a simultaneous fit of v signals and baekgrounds in R? and angular 
eorrelation to the Sun. All these results are eonsistent. 

3.5 Signal Extraction 

In order to test the null hypothesis, the assumption that there are only eleetron neu¬ 
trinos in the solar neutrino flux, the extended maximum likelihood method is used in 
extraeting the CC, ES and neutron (i.e. NC-i-baekground) eontributions in the eandidate 
data set. Baekground eontributions are eonstrained to the measured values diseussed 
above. The undistorted speetrum from Ortiz et a/.il is assumed in the signal de- 
eomposition. Data distributions in Teg, the volume-weighted radial variable {R/RavY 
where i?Av = 600 em is the radius of the aerylie vessel, and cos 6q are simultane¬ 
ously fitted to the probability density funetions (PDFs) generated from Monte Carlo 
simulations, cos 6q is the angle between the reeonstrueted direetion of the event and 
the instantaneous direetion from the Sun to the Earth. This distribution is shown in 
Figure |^a) for the analysis threshold of Teg> 5 MeV and fidueial volume selee- 
tion of i? < 550 em, where R is the reeonstrueted event radius. The forward peak 
(cos 6 ^ 0 ^!) arises from the strong direetionality in the ES reaetion. The cos Oq distribu¬ 
tion for the CC reaetion, before aeeounting for the deteetor response, is expeeted to be 
(l-O. 34 Ocos 6 ' 0 ).ii Fig. P^b) shows the distribution of events in the volume-weighted 
radial variable (R/Ray)^, and Figure p^e) shows the kinetic energy spectrum of the 
selected events. 











H 2 O PDF from the acrylic sources and the radon data 






-'■‘I'- 

A Radon data 

0.05 


H^OPDF 




0.04 

: 

lii 



\ 

0.03 






0.02 






0.01 


^4 


i 



4'' 

A. 

0 




0.5 1 

1.5 2 2.5 i 



(cubic AV radii) 


Figure 13: Comparison between the PDF obtained from the radon spike in the deteetor 
and the H 2 O PDF derived from aerylie souree data. 


The extraetion yields igbT.TlgJI CC events, 263.6t||g ES events, and 
NC events Q The uneertainties given above are statistieal only, and the systematie un- 
eertainties are summarized in Table ^ The extraeted eounts for eaeh neutrino detention 
ehannel above ean be eonverted into integrated fluxes above the kinetie energy thresh¬ 
old of 5 MeV. Assuming an undistorted neutrino speetrumil and using only eross 
seetions for the measured neutrino fluxes from eaeh of the ehannels are (in units of 
10® cm“^s“^): 


0CC = 1.76lo;o5(stat.)l|5;[5^ (syst.) 
0ES = 2.39lo;23(stat.)l|5;j2 (syst.) 
0NC = 5.09+|];43(stat.)+°;4® (syst.). 


The CC and ES results reported here are eonsistent with our earlier results^ for Tes >6.75 
MeV. The exeess of the NC flux over the CC and ES fluxes implies neutrino flavor 
transformations. 


* We note that this rate of neutron events also leads to a lower bound on the proton lifetime for “invisible” 
modes (based on the free neutron that would be left in deuteriumo) in excess of 10^® 


mately 3 orders of magnitude more restrictive than previous limits.0 
mechanism to the solar neutrino NC background is ignored. 


years, approxi- 
The possible contribution of this 








A direct test of the null hypothesis, i.e. there are only in the solar neutrino flux, 
can be readily performed with a simple change of variables: 

0CC ~ 

0NC “f 0/iT 

0ES = 0e + 

where 0e is the flux of the electron component, 0^,- is the flux of non-electron com¬ 
ponent, and 6=0.1559 is the ratio of the elastic scattering cross sections for and v^. 
above the kinetic threshold of 5 MeV. Assuming an undistorted ®B energy spectrum, 
a maximum likelihood extraction using these transformed variables gives the fluxes of 
the electron and non-electron components as: 

(/)e = 1.76l[j;E]^(stat.)l°;°^ (syst.) 

(t)^r = 3.4ll|];4^(stat.)l|5;4^ (syst.) 

Combining the statistical and systematic uncertainties in quadrature, 0^,- is S.dltoJI, 
which is 5.3(7 above zero, and provides strong evidence for flavor transformation con¬ 
sistent with neutrino oscillations.ii@ Adding the Super-Kamiokande ES measurement 
of the ®B fluxi (j)^ = 2.32 ± 0.03(stat.)^Q q® (syst.) as an additional constraint, we find 
0/ir = 3.45lo;62, which is 5.5a above zero. Figure |T^ shows the flux of non-electron 
flavor active neutrinos vs. the flux of electron neutrinos deduced from the SNO data. 
The three bands represent the one standard deviation measurements of the CC, ES, 
and NC rates. The error ellipses represent the 68%, 95%, and 99% joint probability 
contours for 0e and 0^.^- 

Removing the constraint that the solar neutrino energy spectrum is undistorted, the 
signal decomposition is repeated using only the cos 0q and (R/Rav)^ information. The 
total flux of active ®B neutrinos measured with the NC reaction is 

0NC = 6.420}|^(stat.)lo;58 (syst.) 

which is in agreement with the shape constrained value above and with the standard 
solar model predictioni for ®B i/’s, 0 ssm = S.OSlgii- 

3.6 Day-Night Analysis 

The solar neutrino flux results presented above provide strong evidence for neutrino fla¬ 
vor transformation independent of solar model predictions. One of the possible mech¬ 
anism for this flavor transformation is mass-induced neutrino oscillations. For some 



Figure 14: (a) Distribution of cos 6q for R < 550 cm. (b) Distribution of the volume 
weighted radial variable (R/Ray)^- (c) Kinetic energy for R < 550 cm. Also shown 
are the Monte Carlo predictions for CC, ES and NC + background neutron events scaled 
to the fit results, and the calculated spectrum of Cherenkov background (Background) 
events. The dashed lines represent the summed components, and the bands show ±lcr 
uncertainties. All distributions are for events with Teg>5 MeV. 
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Figure 15: Flux of ®B solar neutrinos whieh are /r or r flavor vj'. flux of eleetron 
neutrinos dedueed from the three neutrino reaetions in SNO. The diagonal bands show 
the total flux as predieted by the SSmI (dashed lines) and that measured with the 
NC reaetion in SNO (solid band). The intereepts of these bands with the axes represent 
the ±l(j errors. The bands interseet at the fit values for and indieating that the 
eombined flux results are eonsistent with neutrino flavor transformation assuming no 
distortion in the ®B neutrino energy speetrum. 









Source 

CC Uncert. 

(percent) 

NC Uncert. 

(percent) 

(j)^r Uncert. 
(percent) 

Energy scale f 

-4.2,+4.3 

-6.2,+6.1 

-10.4,+10.3 

Energy resolution f 

-0.9,+0.0 

-0.0,+4.4 

-0.0,+6.8 

Energy non-linearity f 

±0.1 

±0.4 

±0.6 

Vertex resolution f 

±0.0 

±0.1 

±0.2 

Vertex accuracy 

-2.8,+2.9 

±1.8 

±1.4 

Angular resolution 

-0.2,+0.2 

-0.3,+0.3 

-0.3,+0.3 

Internal source pd f 

±0.0 

-1.5,+1.6 

-2.0,+2.2 

External source pd 

±0.1 

-1.0,+1.0 

±1.4 

D 2 O Cherenkov f 

-0.1,+0.2 

-2.6,+1.2 

-3.7,+1.7 

H 2 O Cherenkov 

±0.0 

-0.2,+0.4 

-0.2,+0.6 

AV Cherenkov 

±0.0 

-0.2,+0.2 

-0.3,+0.3 

PMT Cherenkov f 

±0.1 

-2.1,+1.6 

-3.0,+2.2 

Neutron capture 

±0.0 

-4.0,+3.6 

-5.8,+5.2 

Cut acceptance 

-0.2,+0.4 

-0.2,+0.4 

-0.2,+0.4 

Experimental uncertain!) 

' -5.2,+5.2 

-8.5,+9.1 

-13.2,+14.1 

Cross section0 

±1.8 

±1.3 

±1.4 


Table 7: Systematic uncertainties on fluxes. The experimental uncertainty for ES (not 
shown) is -4.8,+5.0 percent, f denotes CC vs. NC anti-correlation. 







values of the mixing parameters, speetral distortions and a measurable dependenee on 
solar zenith angle are expeetedSHS This solar zenith angle dependenee might be 
caused by interaction with the matter of the Earth (the MSW effect) and would depend 
not only on oscillation parameters and neutrino energy, but also on the neutrino path 
length and e~ density through the Earth. 

To look for this distinctive signature of neutrino oscillation, a solar neutrino flux 
analysis similar to that described above was performed for two solar zenith angle 6z 
bins: > 0 (“day”) and 6^^ < 0 (“night”). During the development of this analysis, 

the data were partitioned into two sets of approximately equal livetime (split at July 1, 
2000), each having substantial day and night components. Analysis procedures were 
refined during the analysis of Set 1 and fixed before Set 2 was analyzed with the same 
procedures. The latter thus served as an unbiased test. Unless otherwise stated, the 
analysis presented in the following is for the combined data set. 

Eor each neutrino interaction channel, the asymmetry ratio (^) of the measured day 
flux ((j)D) and night flux 

A, = 2 ~ z = CC,NC,ES 

(PN + (PD 

was determined. Eigure shows the day and night energy spectrum for all the accepted 
events above the kinetic energy threshold of 5 MeV and inside the fiducial volume of 
the inner 550 cm. In the extraction of the neutrino fluxes, backgrounds were subtracted 
separately for the two zenith angle bins. The results were then corrected for the orbital 
eccentricity by normalizing to an Earth-Sun distance of 1 AU. Table || is a summary of 
these extractions, where the day and night fluxes are given for the combined data. A 
consistency test of the six measured fluxes between Sets 1 and 2 yielded a chance 
probability of 8%. A similar test done directly on the three asymmetry parameters 
gave a chance probability of 2%. No systematic has been identified, in either signal or 
background regions, that would suggest that the differences between Set 1 and Set 2 
are other than a statistical fluctuation. Eor the combined analysis, Acc is +2.2a from 
zero, while Aes and Anc are —O.Qcx and —1.2a from zero, respectively. Note that 
Acc and Anc are strongly statistically anti-correlated {p = —0.518), while Acc and 
^Es ip = —0.161) and Aes and Anc (p = —0.106) are moderately anti-correlated. 
Table gives the systematic uncertainties on the asymmetry parameters. 

Systematic checks were made on a set of signals that are continuously present in 
the detector in order to look for any diurnal variation in the detector response. These 
studies include livetime verification using the detector diagnostic triggers (pulsed at 
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Figure 16: (a) Energy spectra (signals -i- background) for day and night bins. The final 
energy bin extends from 13.0 to 20.0 MeV. (b) Difference between the spectra (night 
- day). The day rate was 9.23 ± 0.27 events/day, and the night rate was 9.79 ± 0.24 
events/day. 







Combined 





^(%) 

signal 



4>d 





4>n 





CC 

1.62 

± 

0.08 

± 

0.08 

1.87 

± 

0.07 

± 

0.10 

+14.0 ± 

6.3l|i 

ES 

2.64 

± 

0.37 

± 

0.12 

2.22 

± 

0.30 

± 

0.12 

-17.4 + 

19.5l|^ 

NC 

5.69 

± 

0.66 

± 

0.44 

4.63 

± 

0.57 

± 

0.44 

-20.4 + 

16.911:^ 


Table 8: The results of signal extraction for the two zenith angle bins. The fluxes 
shown here have been normalized to an Earth-Sun distance of 1 AU. An undistorted 
spectrum was assumed in the signal decomposition. The systematic uncertainties 
(combined set) include a component that cancels in the formation of the.4,. Except for 
the dimensionless A, the units are 10® cm“^ s“^. Elux values have been rounded, but 
the asymmetries were calculated with full precision. 




Systematic 

SAcc 

SAes 

SAnc 

Long-term energy scale drift 


■a 

0.2 

Diurnal energy scale variation 

|H 

mm 

1.6 

Directional energy scale var. 

0.2 

1.4 

0.3 

Diurnal energy resolution var. 

0.1 

0.1 

0.3 

Directional energy resolution var. 

0.0 

0.1 

0.0 

Diurnal vertex shift var. 

0.5 

0.6 

0.7 

Directional vertex shift var. 

0.0 

1.1 

0.1 

Diurnal vertex resolution var. 

0.2 

0.7 

0.5 

Directional angular recon. var. 

0.0 

0.1 

0.1 

PMT jJ-'y background 

0.0 

0.2 

0.5 

AV- 1 -H 2 O jS-'y background. 

0.0 

0.6 

0.2 

D 2 O jS-'y, neutrons background. 

0.1 

0.4 

1.2 

External neutrons background. 

0.0 

0.2 

0.4 

Cut acceptance 

0.5 

0.5 

0.5 

Total 

1.5 

2.4 

2.4 


Table 9: Effect of systematic uncertainties on ^ (%). For presentation, uncertainties 
have been symmetrized and rounded. 









Figure 17: Joint probability contours for Atot and Ae- The points indicate the results 
when Atot is allowed to float and when it is constrained to zero. The diagonal band 
indicates the 68% joint contour for the Super-K Aes measurement {Aes = 3.3% ± 
2.2% (stat.)^J' 2 % (sys.).i ) 

5 Hz), variation in the observed muon rate, variation in the detector response to the 
muon-induced secondary neutrons, variation in the detector response to the radioactive 
hotspot on the acrylic vessel, and a neutrino signal extraction based on an east-west 
division instead of a zenith angle division. These studies do not show any significant 
diurnal variation in the detector response. In the study of the diurnal variation of the 
radioactive hotspot, a limit of 3.5% on the rate asymmetry was determined. Because of 
its steeply falling energy spectrum, a 0.3% limit was set for the diurnal variation of the 
detector’s energy scale. 

The asymmetry ratio in the Ve flux Ae and the total active neutrino flux Atotai can be 
readily determined by a change of variables in the signal extraction (see Sec. ^3| ). Fig¬ 
ure 0 shows the Ae vs. Atot joint probability contours. Forcing Atot = 0, as predicted 
by active-only models, yielded the result of Ae = 7.0% ± 4.9% (stat.)tj; 2 % (sys.). 










Region 

XLn/dof 

4>b 

A(%) 

Am? 

tan^ 9 

c.l.(%) 

LMA 

57.0/72 

5.86 

6.4 

5.0 X 10-® 

0.34 

— 

LOW 

67.7/72 

4.95 

5.9 

1.3 X 10-^ 

0.55 

99.5 


Table 10: Best fit points in the MSW plane for global MSW analysis using all solar 
neutrino data. 0 b is the best-fit flux for each point, and has units of 10® cm“^ s“^. 
Am^ has units of eV^. Ae is the predicted asymmetry for each point. 

4 Analysis of Neutrino Mixing Parameters 

Using the day and night energy spectrum in Figure an analysis to constrain the al¬ 
lowed MSW mixing parameters under the 2-z/ flavor mixing framework was performed. 
The radial distribution R and the solar angular correlation (cos6*q) information were 
omitted. The expected number of counts were calculated for each of the energy bins 
by convoluting the ®B neutrino spectrum, il the MSW survival probability, the neutrino 
interaction cross sections and the energy response of the SNO detector. With the high 
energy hep neutrino flux fixed at 9.3 x 10^ cm“^ s“^,l the total ®B flux (ps, the dif¬ 
ference Am? between the squared masses of the two neutrino mass eigenstates, and 
the mixing angle 9 are the only free parameters in a fit of the MSW model to the 
observed data. Figure shows the allowed regions in the Am?-6 space at the 90%, 
95%, 99% and 99.73% confidence levels. 

Additional information from other solar neutrino experiments can be used in con¬ 
structing the allowed Am?-9 regions. In Figure [1^, the solar neutrino fluxes measured 
by the Cl experiment0 and the Ga experiments ,§§11 along with the day and night neu¬ 
trino energy spectra from the Super-Kamiokande experimenti are combined with the 
day and night energy spectra from SNO in a fit of the three free parameters (0 b > Am^ 
and 9). The pp, pep, ^Be and hep neutrino fluxes were fixed at their predicted values 
in the Standard Solar ModelTable summarizes the best fit points in the Am^-9 
plane. It is clear from this table and Figure that this global analysis strongly favors 
the Large Mixing Angle (LMA) region, and excludes the region of tan 9 >\. 

5 Conclusions 

The results that are presented in this paper are truly groundbreaking. The neutral- 
current measurement is the first direct measurement of the total active ®B neutrino flux. 





Figure 18: Allowed regions of the MSW plane determined by a fit to SNO day and 
night energy spectra. The star indicates the best fit. 

The neutrino fluxes measured by this channel provide strong evidence for neutrino fla¬ 
vor transformation, thereby solving the long-standing Solar Neutrino Problem. The 
significance of z/g transforming into or Ur is at 5.3cr level. This neutral-current mea¬ 
surement also verified the Standard Solar Model prediction of the ®B neutrino flux. 

The analysis presented in this paper is also the first direct measurement of the day- 
night asymmetries in the z/g flux and the total neutrino flux. When combining the day 
and night energy spectra from SNO with results from other solar neutrino experiments 
in a global 2-u flavor analysis of the MSW oscillation parameters, the LMA solution is 
strongly favored and the “dark side” (tan 6* >1) is excluded. 

In the coming years, SNO will continue to make significant contributions to solar 
neutrino physics by refining the neutral-current measurement and the day-night asym¬ 
metry with different experimental techniques. It will also attempt to look for other 
possible signatures of neutrino oscillation (e.g. spectral distortion) and other rare pro¬ 
cesses (e.g. nucleon decays, anti-z/) in its future physics program. 





Figure 19: Allowed regions of the MSW plane determined by a fit to SNO day 
and night energy speetra, the solar neutrino fluxes from the Cl and Ga experiments and 
the day and night energy speetra from the Super-Kamiokande experiment. The star 
indieates the best fit. 
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